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We calculate the spectral properties of the one-dimensional Holstein and breathing polarons using 
the self-consistent Born approximation. The Holstein model electron-phonon coupling is momentum 
independent while the breathing coupling increases monotonically with the phonon momentum. We 
find that for a linear or tight binding electron dispersion: i) for the same value of the dimensionless 
coupling the quasiparticle renormalization at small momentum in the breathing polaron is much 
smaller, ii) the quasiparticle renormalization at small momentum in the breathing polaron increases 
with phonon frequency unlike in the Holstein model where it decreases, iii) in the Holstein model the 
quasiparticle dispersion displays a kink and a small gap at an excitation energy equal to the phonon 
frequency ljo while in the breathing model it displays two gaps, one at excitation energy uJo and an- 
other one at 2ijjq. These differences have two reasons: first, the momentum of the relevant scattered 
phonons increases with increasing polaron momentum and second, the breathing bare coupling is an 
increasing function of the phonon momentum. These result in an effective electron-phonon coupling 
for the breathing model which is an increasing function of the total polaron momentum, such that 
the small momentum polaron is in the weak coupling regime while the large momentum one is in 
the strong coupling regime. However the first reason does not hold if the free electron dispersion has 
low energy states separated by large momentum, as in a higher dimensional system for example, in 
which situation the difference between the two models becomes less significant. 



I. INTRODUCTION 

Interest in the spectral properties of strongly coupled 
electron-phonon systems has increased due to the discov- 
ery of a kink in the quasiparticle dispersio niiMi'^'^i mea- 
sured by angle-resolved photoemission (ARPES) exper- 
iments in cuprate superconductors. The energy of this 
kink is characteristic of a specific optical phonon mode. 
Furthermore, a peculiar isotope-effect, manifesting in a 
shift of the quasiparticle dispersion at high energy upon 
oxygen isotope substitution, was recently measured in 
the high Tc materials^, undoubtedly showing that the 
electron-phonon interaction plays an important role in 
the spectral properties of cuprates. 

The electron-phonon interaction has been investigated 
in a variety of model o^i^i^ . Although the ground state 
properties of these models received great attention, less 
was paid to the spectral features. Noticeable excep- 
tions include the HolsteinSiiS and the Frolichii models 
which are characterized by a momentum independent or a 
weakly momentum dependent electron-phonon coupling. 
However, in many systems, such as cuprates supercon- 
ductors or organic materials, the electron-phonon cou- 
pling is strongly momentum dependent. This can give 
rise to distinctive properties of the single particle spectral 
features, as we will discuss in this paper. In particular we 
will consider a sinusoidal dependence of the coupling on 
the phonon momentum, which is realized in many sys- 
tems including the SSH coupling in polyacetilens^ and 
the coupling to the breathing mode in cuprate supercon- 
ductorsiajaii. 

The treatment of a multi-electron system coupled with 
phonons is an extremely complex and difficult problem. 



In general, phonons mediate an effective attraction be- 
tween electrons and consequently the system becomes 
susceptible to various kinds of instabilitie3i^*i&, which 
might have a significant effect on the photoemission spec- 
tra. However, in this paper we investigate the effect of 
phonons on a single quasiparticle only, thus neglecting 
electron-electron or electron-hole scattering processes at 
the Fermi surface. We believe this to be suitable in de- 
scribing the properties at energy or temperature scales 
larger than those associated with the ground state insta- 
bility, e.g. the normal state of superconductors. 

In fact, as different investigations have showniStii, the 
polaron models, which consist of a single electron inter- 
acting with phonons capture much of the physics seen 
in the ARPES experiments on materials with significant 
electron-phonon interactionSiiS. This includes the kink 
in the quasiparticle dispersion observed at the phonon 
characteristic frequency. One expects that the single 
particle description of the influence of electron-phonon 
coupling would be valid to describe the photoelectron 
or inverse photoelectron spectral function in insulators 
and semiconductors in which the bands are either full 
or empty. However things are more complicated for a 
strongly correlated insulator in which the conductivity 
gap is a result of electron correlations and the material 
is an insulator in spite of having a half filled band. In 
systems described by the Hubbard or t-J models the sin- 
gle hole spectral function will be infiuenced by the in- 
teraction with spin fluctuations. This results in strongly 
dressed quasiparticles even without the electron-phonon 
coupling, as calculations employing self consistent Born 
approximation (SCBA)^" and exact diagrammatic quan- 
tum Monte Carlo (QMC)2i shows. The later calculation 
also shows that the SCBA, which neglects all the magnon 
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crossing diagrams, is a good approximation to the quasi- 
particle mass renormalization and its dispersion. When 
the phonons are considered, in the simplest approxima- 
tion we can neglect all the crossing diagrams involving 
phonons and magnons and treat the interaction with 
phonons in a single particle fashion by taking the hole 
band dispersion as that given by the t-J model calcula- 
tions. The approximations involving the polaron models 
under investigation here are done in this spirit. 

In this paper we investigate the spectral properties of 
two polaron models: The one-dimensional Holstein (H) 
model® and a one-dimensional version of the breathing 
model (B)i£, relevant for cuprates supcrconductors22i2^. 
The electron-phonon coupling in the H Hamiltonian is a 
constant, thus independent of the phonon momentum. 
The B model has a sinusoidal momentum dependent 
electron-phonon coupling, being small (large) for scat- 
tering with small (large) momentum phonons. In both 
cases we consider dispersionless optical phonons with fre- 
quency loq. Although for the sake of simplicity our cal- 
culations are done in one dimension, the conclusions and 
the qualitative properties of the spectra are independent 
of dimensionality. 

The difference between the properties of the two mod- 
els emphasizes the importance of momentum dependent 
couplings. The specific momentum dependence of the B 
bare electron-phonon coupling results in an effective cou- 
pling which is an increasing function of the total polaron 
momentum, such that the small momentum polaron is in 
the weak coupling regime while the large momentum one 
is in the strong coupling regime. This might be germane 
to the peculiar behavior of the high energy quasiparticle 
dispersion^ or the temperature dependency of the pho- 
toemission linewidth24 in cuprate superconductors. 

The method we use in our calculations is the self- 
consistent Born approximation (SCBA)'^^. Although this 
method is an uncontrolled approximation which neglects 
electron-phonon vertex corrections in the self-energy cal- 
culation, the results at small coupling are in good agree- 
ment with numerically exact quantum Monte Carlo cal- 
culations. Because the contribution of configurations 
with an infinite number of phonons is considered, this 
approximation manages to capture the physics at large 
momentum, such as the gap and the flattening of the 
quasiparticle dispersion at the phonon energy^^. It can 
easily provide information about excited states and thus 
determine the polaron spectral properties, unlike most 
QMC methodsiii2L2& fit to calculate only the ground 
state properties. 

The paper is organized as follows. In Sec. Ill Al we 
introduce the H and B Hamiltonians. The dimension- 
less electron-phonon coupling is defined in Sec. Ill Bl The 
SCBA method is discussed in Sec. Ill CI The results are 
presented in Sec. lIIII and their significance is discussed in 
Sec. II VI A short summary and the conclusions are given 
in Sec. |V| 



II. FORMALISM 
A. Model 

The Holstein Hamiltonian in real space is given by 
Hh = X! {cUj + h.c.) + uJo^blbi+g^n.iXi , (1) 

<ij> i i 

where Ci (bi) is the electron (phonon) annihilation opera- 
tor at site i. The first term describes the electron kinetic 
energy. The second term describes a set of independent 
oscillators with frequency ujq at every site. The electron- 
phonon coupling in the Holstein model is local and is 
described by the last term of Eq. 2] where the electron 
density Ui = c|ci couples with the lattice displacement 

= vA^^^l + ^i*^ strength g. 

Our breathing Hamiltonian is a one-dimensional ver- 
sion of the model which describes the coupling of Zhang- 
Rice (ZR) singlcts^^ with the Cu-0 bond-stretching vi- 
brations in high Tc superconductor o^^i^? . In cuprates a 
ZR singlet is a bound state between a hole on the Cu and 
a hole on the four neighboring O atoms. It's energy is sta- 
bilized by the Cu-0 hybridization term and therefore is 
influenced by the Cu-0 distance. In our one dimensional 
model we consider a set of independent, in-between sites 
oscillators (the analogue of the O atoms) which modu- 
late the charge carrier's (the analogue of the ZR singlet) 
on-site energy. Therefore we define the B Hamiltonian as 

Hb = -tE<„->(clc, + h.c.) + u;o bl^K+i (2) 

Both the H and the B model can be written in the mo- 
mentum representation as 

HiH,B) = <^)4cfc + J2 + (3) 

k q 

■4^ J2 ^(h,b) {q)cl_gCk{bl + b^q) 

where 

is the H electron-phonon coupling and 

is the B electron-phonon coupling. Notice that in the mo- 
mentum representation the H coupling is a constant and 
the B coupling is an increasing function of the phonon 
momentum for small momenta. 

In Eq.^and Eq.|21the free electron part of the Hamilto- 
nian was introduced as a tight-binding hopping. However 
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in order to study the influence of electron dispersion on 
the polaron properties in this paper we also employ cal- 
culations which consider different forms of electron dis- 
persion. 



B. Electron-phonon coupling 

In this paper, we define the dimensionless electron- 
phonon coupling as the ratio of the lattice deformation 
energy of a localized electron and the kinetic (delocaliza- 
tion) energy of the electron. This definition is encoun- 
tered mostly in polaron studies focusing on the transition 
from large to small polarongSSi2£. For the Holstein model 
the deformation energy is given by (see Eq. HTH)! 



E = : 2— 

^" 2MuJo LOO ' 
and the dimensionless coupling is defined as 
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where the half-bandwidth W/2 is taken as a measure of 
the free electron kinetic energy. For a simple tight bind- 
ing dispersion W — 2zt, where z is the coordination num- 
ber. 

Analogously, for the breathing model the lattice defor- 
mation energy is (see Appendix I A 2|l 
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and the dimensionless coupling is hence defined as 
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Compared to the Holstein case, the lattice deformation 
energy for the breathing model has an extra factor of z 
as can be seen in Eq. |S1 which results form the fact that 
the electron interacts with z neighboring oscillators. 

There are several other definitions for the dimension- 
less electron-phonon coupling throughout the literature. 
For instance in the BCS theory the definition of the di- 
mensionless coupling is A = VN{0), where V and A^(0) 
are the effective electron-electron attraction and respec- 
tively the density of states at the Fermi level. In the 
weak coupling regime the effective interaction at small 
frequency or in the antiadiabatic limit {ujo/t oo) is 
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If the integration over all q momenta is considered for 
the breathing model, then 
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When choosing the density of states A'^(O) = 1/W, this 
definition of the dimensionless coupling coincides with 
our definition. 

In the Migdal-Eliashberg theory of superconductivity, 
)^ME _ 2a'^F{ujo)/uJo where only scattering processes at 
the Fermi surface are considere d?^'?^ in a^Floj). This 
definition is equivalent tc^^ 
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where to* is the quasiparticle renormalized effective mass 
at the Fermi surface. A*^^ can be directly determined 
from experiments since a'^F{Lu) and to * /to can be mea- 
sured in tunnelingSi, neutron^ or respectively in the 
ARPES experiments^. In our case, the effective mass at 
the bottom of the band (which is the zero energy state 
of our system and in many respects similar to the Fermi 
surface in a many-electron system) can be determined 
from the equation 

{E{k)-Em{l~ — {0,E{0)) = e{k)-e{0) + — {0,Emk 

(14) 

valid at small k. The ratio of the effective mass and the 
un-renormalized mass is 



TO _ E{k) - E{0) 
m* e{k) — e(0) 

Therefore, we have 
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This definition directly relates the quasiparticle weight 
Zq to the coupling constant A*^^. In the first order per- 
turbation theory one gets 
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By inspecting Eg. llTl one can see that the most important 
contribution to the integral comes form small q where the 
denominator is small. For the B model this will introduce 
a significant difference between and our Ab (Eq. O . 
Unlike our definition which assumes average over all pos- 
sible phonon momenta, A;^^^ is determined by small q 
phonon scatterings which are characterized by small 7(g) 
in the B model. 

The definition of A*^^ used in the Migdal-Eliashberg 
theory was intended to describe the properties at the 
Fermi level or at small energy excitation. It does not 
properly describe the physics at larger energy since the 
most relevant scattering processes in this case imply 
larger phonon momenta as we will discuss in the next 
two sections. 



C. Method 

In the SCBA the electron self-energy is obtained by 
summing over all non-crossing diagrams. In this approx- 
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FIG. 1: Schematic representation of the Dyson equation (up- 
per). In the SCBA approximation the self-energy is a sum- 
mation of aU non-crossing diagrams (middle) and can be de- 
termined self-consistently using Eg. \'20\ ( bottom). 



imation the calculation of the self-energy can be reduced 
to the following set of equations: 



G(k,u) 



oj — e(fc) — I](fc, Lo) ' 



D(q,^l) 



(18) 



(19) 



and 



Y.(k,uj) 



(27r)2 



dn J dq \j{q)\^D(q,n)G{k-q,oj^n) , 

(20) 

where G{k,uj) and D(q,i}) are the free electron and re- 
spectively the free phonon Green's functions (see also 
Fig. . The frequency integration in Eq. [201 can be ex- 
plicitly completed and one gets 



^{k,i^) = [ dq- 
2n J L 



luq — e(fc ~ <?) ~ S(fc — q^uj — ujq) 

(21) 



where is the phonon frequency. Eq. [2] is solved self- 
consistently and the spectral function is determined as 
A(k,uj) — — ilmG(fc,Li;). The overlap of the polaron 
state \vk > (the lowest energy state at a particular mo- 
mentum k which in general is an isolated pole in the 
Green's function) with the free electron state, c^|0) 
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is called the quasiparticle weight and can be calculated 
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FIG. 2: Numerically exact QMC (diamonds) and SCBA 
(squares) results for polaron energy Eo (a) and quasiparti- 
cle weight Zo (b) versus the dimensionless electron-phonon 
coupling A (see Eq. at = 0. Holstein model with t = 1 
and LOO = 0.1 t and 2M = 1. 



The neglect of vertex corrections in the SCBA results 
in a failure of this approximation in the intermediate and 
strong coupling regime. This is apparent in Fig. |2] where 
a comparison of SCBA and exact QMC results for the H 
model energy and quasiparticle weight at A: = is shown. 
In fact the polaron properties at fc = are better approx- 
imated in a regular Rayleigh-Schrodinger (RS) perturba- 
tion theory than in SCBA'^^. However, at large k where 
the polaron energy approaches the phonon frequency, the 
RS perturbation theory fails and the SCBA which con- 
siders configurations involving a large number of phonons 
provides good results. In order to show that, we compare 
the SCBA results with the exact ones obtained by using 
the diagrammatic quantum Monte Carlo technique^SiSi 
in Fig.|3| As can be seen in Fig.|21-a, the dispersion be- 
comes flat when the polaron energy gets close to the the 
phonon frequency Eq+luo: Eq being the bottom of the po- 
laronic band. The flattening is accompanied by a strong 
reduction of the quasiparticle weight, as shown in Fig.|2| 
-b. The SCBA slightly underestimates the energy correc- 
tions in the flat dispersion region but otherwise captures 
all these features successfully. 

The H polaron physics at small coupling has been in- 
vestigated with different techniques for many years and 
is rather well understood now. The polaron character 
changes from almost free electron at the bottom of the 
band to a one phonon plus one electron state in the flat 
dispersion region where almost all the momentum is car- 
ried by the phonon. The states below the one phonon 
threshold Eq + ujq are bound statesii, characterized by 
delta peaks in the spectral representation. Above it, 
there is a continuum of states and the electron self-energy 
acquires a finite imaginary part. However, at large en- 
ergies relative to the self-energy's imaginary part is 
small and the spectral representation is characterized by 
sharp peaks with a dispersion close to the free electron 
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FIG. 3: Holstein polaron dispersion E{k) (a) and quasi- 
particle weight Z{k) (b) calculated with SCBA (dotted line) 
and QMC (circles), c) Spectral representation obtained with 
SCBA. Everywhere \h = 0.25 and ujq = 0.5 t. 



one. In a photoemission experiment the dispersion of the 
high intensity peak would exhibit something resembUng 
a kink, although a small gap at the phonon energy and 
sudden onset of broadening appears at the same time. 
The spectral function at momenta corresponding to en- 
ergies well below the phonon frequency will consist of a 
sharp peak followed by a broad satellite, such as it is ob- 
served in surface studies of Be^^ for instance. This can 
be seen in Fig. 13 c) where a false color plot of spectral 
intensity obtained with SCBA is shown (see also Fig. [S] 
-a). The physics discussed in this paragraph (i.e. the dis- 
persion flattening and the gap in the polaron spectrum) 
is presumably true even for very small electron-phonon 
couplings. However, for such small couplings these fea- 
tures are too small to be captured experimentally or even 
numerically. 
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FIG. 4: H (circles) and B (diamonds) polaron energy (a) and 
quasiparticle weight (b) versus A at fc = and too = 0-2t. The 
open symbols are the results obtained after only one iteration 
(first order in D). The dashed line in (b) corresponds to 
^""^ = - 1 (see Eq.CSll. 
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III. RESULTS 

In our calculations we take t — 1 and 2M = 1. For the 
one-dimensional case the coordination number z = 2. 

In Fig. ^ -a and -b we show the polaron energy Eq 
and respectively the quasiparticle weight Zq versus the 
dimensionless coupling constant A at fc = 0. For the H 
model, the decrease of Zq with increasing A is not very 
different from the one given by the ME theory (dashed 
line in Fig.0[b given by Ea. [TC|) . showing that the ME def- 
inition of A is similar to ours. However, for the B model 
the quasiparticle weight Zq and the polaron energy Eq 
decrease much slower with increasing A. As mentioned 
in the previous section this is due to the small momentum 
scattering processes implied in the renormalization of Zq 
(see Eq. I17|) and, similarly, in the determination of the 
self-energy. It is worth pointing out that, unlike the H 
polaron, even for values of A « 1, the B polaron remains 
in the weak coupling regime and hence the difference be- 
tween the fully convergent SCBA (full symbols) and the 



FIG. 5: The quasiparticle weight Zo (a) and the polaron en- 
ergy Eo (b) at zero momentum versus phonon frequency uio for 
the H and B models at \h = 0.2 and respectively As = 1.3. 



first order perturbation theory (i.e. only the first SCBA 
iteration, empty symbols) is small. Another interesting 
feature is that for the same value of Zq the B polaron 
energy is lower than the H one, showing that the ratio 
between the energy renormalization and the quasiparticle 
weight renormalization is different for the two models. 

Another important difference between the two models 
is the dependence of the polaron properties at the band 
bottom on the phonon frequency uq. While for the H 
case an increase of loq results in an increase of the quasi- 
particle weight Zq the opposite behavior is seen for the B 
model. This is illustrated in Fig. [SJ-a. The reason for the 
reduction of Zq with increasing ujq in the B model can be 
easily understood by noticing (see Eq. 117(1 that a larger 
value of loq reduces the importance of the q dependence in 
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FIG. 6: Quasiparticle weight Z{k) (a) and the polaron energy 
E(k) (b) for the H and B model for wo = 0.2 t at coupling 
Xh = 0.2 and Xb = 1.3 respectively. The dotted horizontal 
lines in (b) marks the first phonon threshold energy Eo + ujo- 



the polaron properties calculation. As discussed earlier, 
the strong momentum dependent coupling is responsible 
for the weak Zq rcnormalization of the B polaron and 
thus an increase of loq would result in a larger effective 
coupling and implicitly in a smaller Zq. 

Aside from the different A and uq dependency of the 
two models at zero momentum, the momentum depen- 
dent properties also exhibit different behaviors. This is 
shown in Figs. El El and El where the k dependent 
properties for the two models are illustrated. 

In Fig. Eland Fig. da and -b we have chosen the value 
of A such that both models yield the same quasiparticle 
weight at the bottom of the band. Thus, the choice of 
Xh = 0.2 and Xb = 1-3 results in Zoh = Zob = 0.83, im- 
plying that both models are in the weak coupling regime. 
As mentioned earlier, the B polaron energy at fc = is 
lower. At large k, just below the first phonon thresh- 
old energy Eq + loq, both polarons display a flat disper- 
sion and a reduced quasiparticle weight. However, the 
B polaron quasiparticle weight at large k is substantially 
larger (for instance Zsik = it) Ki 0.014 within a numer- 
ical precision of 10"'^ and Znik — n) < 10"^ ), making 
the B polaron state at large k distinguishable in the spec- 
tral plot (Fig. [7|-b) in contrast to the H one (Fig. Qa). 
At energies larger than Eq + ivo the spectral intensity of 
the B quasiparticle is much smaller than the H model 
one, unlike the situation at the band bottom where both 
models have the same Zq. This large momentum behav- 
ior points to a stronger effective coupling for the B model 
at large k. 

While the dispersion of both models displays a gap at 
Eq+loq, the B polaron shows a second gap at the second 
phonons threshold energy Eq + 2ujq. This can be seen 
in Fig. |71b but also occurs for smaller values of the di- 
mensionless coupling as shown in Fig. Q-c for the value 
of Xb = 0.5. This value was chosen such that the ground 
state energy Eq of the B model is equal to that of the Hol- 
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FIG. 7: Spectral representation for; a) H model. Ah = 0.2, 
tight binding dispersion, b) B model, Xb = 1.3, tight binding 
dispersion, c) B model, Xb ~ 0.5, tight binding dispersion, d) 
B model, Ag" = 1.6, linear dispersion with vf ~ t f Eg. 1251 . 
e) H model, A^" ~ 0.9, linear dispersion with vf ~ t, f) 
H and B model. Ah ~ Xb ~ 0.2, with dispersion e{k) = 
— 2tcos(2fc). Everywhere loq — 0.2t. The arrows indicate the 
first (lower) and the second (upper) phonon threshold energies 
corresponding to Eq + ujq and respectively Eq + 2ijjo. 



stein polaron one shown in Fig. da. The situation can be 
even clearer visualized by comparing Fig.Ela with Fig. El 
b, where the energy distribution curves (EDC) for H and 
respectively B cases are shown. 

An even more interesting effect is noticed if a linear 
dispersion for the free electron is considered 



e(k) = VF\k\ , 



(24) 
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FIG. 8: Spectral function (energy distribution curves) for a) 
H model with Xh = 0.2 and b) B model with As = 0.5. Free 
electron tight binding dispersion and ujo — 0.2t is considered. 




FIG. 9; Spectral function (energy distribution curves) for a) 
H model with A^" = 0.9 and b) B model with A^" = 1.6. 
Free electron linear dispersion {vf = t) and loq — 0.2t is 
considered. 

with a value oi vp close to one or larger. In this case one 
can take vp to be a measure of the free electron kinetic 
energji^ and thus define the dimensionless coupUng as 



small coupling strength 7(9), originating from the rapid 
increase of the electron energy with k. At larger energy, 
close to Eq + 2^0, the relevant phonon momenta q im- 
plied in the scattering are larger and the physics is con- 
sequently determined by a larger effective coupling. As 
a result a noticeable kink appears at this energy in the 
spectrum. This effect (i.e. kink at Eq -f 2luo but no no- 
ticeable one at Eq + wo) is a result of an electron-phonon 
coupling which is an increasing function of the polaron 
momentum and it is hence not seen in the Holstein model 
even for the case of a linear electronic dispersion (see 
Fig. EJe and the corresponding EDC plot in Fig. |5[a). 

The differences between the two models discussed 
above are a consequence of two effects: i) the strong q 
dependence of the bare electron-phonon coupling 7(5) in 
the B model and ii) the polaron properties at small k 
are most strongly influenced by the small momentum q 
phonons. However the second statement is not true if the 
free electron dispersion has low energy states separated 
by large q as we will discuss in the next section (Sec lIV|l . 
In order to show this we choose a free electron dispersion 

e(/c) = -2icos(2fc) , (26) 

which is double degenerate with the lowest energy values 
at k = and k — tt respectively. We find that for this 
electronic dispersion the differences between the H and 
the B model are very small, less than 0.1%, and therfore 
not discernible in the spectral representation plot shown 
in Fig.Qf. 

The general features of the polaron spectral function 
illustrated in Figs. |H1 and |3 show a remarkable resem- 
blance with the photocmission data in materials with sig- 
nificant electron-phonon interactions^^. For momenta 
which correspond to energies below the phonon frequency 
one can see a sharp peak followed by a broad satellite. 
At large momenta, the broad satellite found at energies 
well above the phonon frequency sharpens and follows 
the free electron dispersion while the intensity of the 
sharp peak below loq vanishes. We want to point that 
this picture, with the spectral function determined by 
two main branches, is quite different form the one in- 
correctly presented in some popular textbooks^S where 
the spectrum is characterized by a single quasiparticle 
peak which changes its dispersion slope when crossing 
the phonon energy threshold. 



\ lin 



E.. 



PB 



1 



vp 2MujQ vpLOo 



(25) 



The resulting B polaron dispersion is shown in Fig. [3-d 
(see also the corresponding EDC plot in Fig.lHl-b). While 
it displays a gap at Eq + 2cjo, no distinguishable gap or 
kink can be seen at the first phonon threshold energy 
Eq + luq. This free-electron like behavior of the polaron 
at Eq + ujq is due to the fact that the physics there is 
determined by very small q scatterings, characterized by 



IV. DISCUSSION 

For a tight-binding or a linear dispersion, the B polaron 
properties at small momentum are determined by small 
q scatterings characterized by small jiq). Let's consider 
the first order self-energy diagram contribution 



S(fc,w) ^ J dq 



uj — ujQ — e(fc — q) 



(27) 
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The most significant contribution below the one phonon 
threshold, i.e. when to < e(0) + Wqi is given by the values 
of transferred momentum q for which e{k—q) is minimum 
which implies q Ki k. The value of k where the polaron 
energy reaches the one-phonon threshold is small (fci in 
Fig.llOI-a) and thus, the relevant scatterings which deter- 
mine the physics below Eq~\-ujq occur at small momentum 
g, implying a small effective coupling. 

For a linear dispersion with vp ^ 1 ov larger, the value 
of k where the first phonon threshold energy is reached is 
even smaller. Besides, the electronic density of states at 
small momentum, which is proportional to the number 
of relevant scatterings, is also much smaller then for a ID 
tight-binding dispersion^^. These two conditions yield a 
very small effective coupling at the first phonon thresh- 
old. Consequently, even for large values of A as defined 
in Eq. 1251 the polaron dispersion exhibits an extremely 
narrow gap, hardly discernible in our calculation. 

The B model properties at larger momentum indicate 
an increase of the effective electron-phonon coupling with 
increasing the polaron momentum, due to the increase 
of the momentum of the relevant scattered phonons. 
For example, the most significant scattering at the two 
phonon threshold {lo k, Eq + 2wo) are those for which 
e{k — q) K, ujQ. If fci is the value of the electron momen- 
tum for which e(fci) = e(— fci) = wq the resulting values 
for the phonon momentum q relevant in the scattering 
are q fa k ± ki. The larger q solution, i.e. q k + ki {k2 
in Fig. Iim-a). implies a larger effective coupling. 

A similar analysis of the higher order self-energy dia- 
grams leads to the same conclusion: The effective cou- 
pling in the B polaron model is increasing with momen- 
tum, and hence the small energy and momentum B po- 
laron is characterized by a small effective coupling while 
the high energy and large momentum properties are de- 
termined by a large effective coupling. Moreover, at low 
energy the contribution of the crossing diagrams is com- 
parable to the contribution of the crossing diagrams in 
a H model at small coupling. Thus, the SCBA solution 
is a good approximation for the low energy, small k, po- 
laron. Nevertheless, the large momentum properties are 
characterized by large effective couplings and hence the 
SCBA approximation is questionable in that region. 

The conclusions based on the analysis of self-energy 
diagrams discussed above implicitly assume that a small 
e{k — q) implies small q, when the polaron momentum 
k is small. However, as we have shown in the previous 
section, if the electron dispersion has low energy states 
separated by large q the above assumption is invalid. 
Thus the effective coupling at small energy is substan- 
tially larger and besides there is no significant increase of 
the effective coupling with polaron momentum or energy. 

For example the strong dependence of the effective cou- 
pling on the momentum does not hold in a metallic sys- 
tem with a large Fermi surface. To see this, one can think 
of the Fermi surface as a degenerate ground state for the 
quasiparticle^. Unlike the non-degenerate ground state 
case where the small energy scattering is restricted to 



a) 

E(kjK 



one electron +'\ 




FIG. 10: a) Energy levels for zero coupling {g = 0). The solid 
line shows the non-interacting electron dispersion, the dashed 
ones are one phonon + one electron states and the dotted one 
two phonons + one electron, etc. k\ is the momentum where 
the electron dispersion reaches and has the most significant 
contribution to the scattering at the first phonon threshold. 
k2 is the value of momentum with significant contribution to 
the scattering at the second phonon threshold, b) For the po- 
laron with momentum k the scattering with momentum q has 
a significant contribution. Although the polaron momentum 
is close to kp, thus ka = k — kp is small, q = ka + ki, is large 
because kb which connects two points on the Fermi surface 
can be large. 

small momentum, the relevant scatterings in the pres- 
ence of a Fermi surface are restricted to small momen- 
tum plus a vector which connects two points on the Fermi 
surface {q = ka + fcb in Fig.^|b)). Therefore scatterings 
with large momentum, for which 7(g) is large, are rele- 
vant. This makes the effective coupling to be large even 
at small energy. Because now at both low energy and 
high energy there are relevant large momentum scatter- 
ings the dependence of the effective coupling strength on 
the energy and momentum will be modest. 

Although dimensionality plays an important role in de- 
termining the quantitative polaron propertie a^^i"^" , the 
main features of the B polaron are a consequence of 
strong momentum dependence of the bare electron- 
phonon coupling. Therefore the main conclusions of this 
study remain valid for 2D or 3D systems as long as the 
electron-phonon coupling has a similar q dependence (in- 
creases with q). 

The particularities of the B polaron discussed in this 
paper might be relevant for the weakly doped cuprates 
characterized by very small Fermi pockets around the 
(±7r/2, ±7r/2) points in the Brillouin Zone which makes 
the situation very similar to the one captured by the po- 
laron model. The quasiparticles in this system have a 
four-fold degenerate ground state with the states sepa- 
rated by (0, tt) or (tt, 0) vectors^S. The relevant phonon 
mode is believed to be the half-breathing one for which 
the electron-phonon coupling has a strong q dependence, 
being small at small q and large at (0, tt). If for some rea- 
sons the scattering between the (±7r/2, ±7r/2) points is 
restricted or if its importance is small because of an addi- 
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tional k dependence in the bare electron-phonon couphng 
"f{k,q) the formed polarons wiU be characterized by a mo- 
mentum and energy dependent effective couphng. Conse- 
quently, this implies interesting energy and temperature 
dependent properties such as a strongly temperature de- 
pendent quasiparticle photoemission linewidthiS4. How- 
ever, if the scattering between the (±7r/2, ±7r/2) states 
is relevant, the q dependence of the coupling will not be 
very important and presumably a Holstein like coupling 
can be used to describe the physics as well, analogous to 
the ID case discussed in Fig. [7[f. 
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APPENDIX A: EFFECTIVE COUPLING 
1. Holstein 



V. CONCLUSIONS 

In this paper we study the spectral properties of the 
one-dimensional Holstein and breathing polaron mod- 
els. While the H model has a momentum independent 
electron-phonon coupling, the B coupling is a mono- 
tonic increasing function of the phonon momentum. We 
find the renormalization of the quasiparticle properties 
at small momentum in the B model weaker than that in 
the H one for a dimensionless coupling constant defined 
as the ratio of the lattice deformation energy and the 
kinetic energy of the electron. With increasing phonon 
frequency wq, the B quasiparticle weight decreases in con- 
trast to the H case where it increases. The quasiparticle 
dispersion in the H model displays a kink and a small gap 
at an excitation energy equal to luq. On the other hand, 
the B model dispersion, besides the gap at loq, exhibits 
another, more pronounced one at energy 2luo. This is 
due to the following facts: the momentum q of the rel- 
evant phonons in the scattering process increases with 
increasing the total polaron momentum and that the B 
coupling is an increasing function of q. As a consequence 
the renormalization of the B polaron increases with in- 
creasing polaron momentum. However the first fact is 
dependent on the free electron dispersion and does not 
hold if it has low energy states separated by large mo- 
mentum. In that case the difference between the B and 
the H polaron becomes less significant. 

The specific momentum dependence of the B bare 
electron-phonon coupling results in an effective coupling 
which is an increasing function of the total polaron mo- 
mentum, such that the small momentum polaron prop- 
erties are weak coupling while the large momentum ones 
are strong coupling. This might be relevant for explain- 
ing the peculiar behavior of the high energy quasiparticle 
dispersion^ or the temperature dependency of the pho- 
toemission linewidtli24 in cuprate superconductors. 
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The last two terms of the Holstein Hamiltonian 
(Eq.^ can be diagonalized via the Lang-Firsov canonical 
transformation^ defined by a unitary operator e'^ where 



S = 



lU()\/2Mluo 
Using the expansion 



(Al) 



A+[S,A] + -[S, [S,A]] 



(A2) 



we find for the transformed phonon annihilation opera- 
tor: 
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and similarly for the electron operator 



(A3) 



(A4) 



The tilde mark is used to label the transformed operators. 
The Holstein Hamiltonian in the new basis can be written 



as 



H^Ho + Ht 



Ho = uJoblbi - -r-r^^ nf 



2MujI ^ 



Ht = ~tY, i^c,x}X, + H.c. 

<ij> 



(A5) 
(A6) 
(A7) 



where 



(A8) 



For a single electron the density term in Eq. IA6I sim- 
plifies, yielding 



^0 = y^^o^i^t 



2MujI 



(A9) 
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One can define a dimensionless coupling constant for 
the electron-phonon interaction as the ratio between the 
gained lattice deformation energy (see second term in i/o) 



77^ £ 

2MuJo ujo 



(AlO) 



and the bare electron kinetic energy, taken to be the half- 
bandwidth, W/2 = zt, 



and for the breathing-Hamiltonian 



Ho = y^^o6-fc 



m 9_ 



2MujI 



(A17) 



<ij> 



(A18) 



2Ep^ 



1 
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2MuJo zLUot 



(All) where 



2. Breathing 

The electron-phonon interaction term in the B Hamil- 
tonian (see last term in Eq. can, by changing the order 
of summation, be re-write it in the form 



(A19) 



In the case of a system only containing a single electron 
the density terms contained in Eo . I A 1 71 simplify, yielding 



Hint = g'^Xi_i{ni - Ui-i) , 



(A12) 



where Xi is the lattice displacement. The last two parts 
of the Hamiltonian are once again diagonalized via the 
unitary operator e"^ where 



Ho = y^cjQ&j bj 



2g^ 



2MojI 



(A20) 



Hence, the lattice deformation energy is found to be 



S 



ujo\'2A'iLJo 2 2 



(A13) 



Using the expansion of Ea. lA2l we find for the transformed 
phonon and electron annihilation operators: 



3 



Wo 



y/2M~ 



Wo 



— r^=ib^ T~b. 1 ) 



(A14) 



(A15) 



Substituting these transformed operators we find in the 
new basis 



'2 y/2MuJo 

= — (h^ 1 + &, 1 ) H ^-77(71, — n,_i ) 

iMM. 
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(A21) 



In contrast to the H model (Eq. lAlOp we incur an addi- 
tional factor of 2 due to a coupling with two neighbor- 
ing oscillators. For higher dimensionality the number of 
neighboring oscillators is given by the coordination num- 
ber z, thus replacing the factor 2 in Ea. IA2l| with z. 

We thus define the dimensionless coupling constant for 
the electron-brcathing-phonon interaction 



2B 



1 



2Mloo zcuot 



(A22) 



This definition is identical to that found for the Holstein 
polaron besides the overall factor of z. 
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